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Abstract: Peer to peer (P2P) systems are moving from appHcation specific architectures to a generic 
service oriented design philosophy. This raises interesting problems in connection with providing 
useful P2P middleware services that are capable of dealing with resource assignment and manage- 
ment in a large-scale, heterogeneous and unreliable environment. One such service, the slicing ser- 
vice, has been proposed to allow for an automatic partitioning of P2P networks into groups (slices) 
that represent a controllable amount of some resource and that are also relatively homogeneous with 
respect to that resource, in the face of churn and other failures. In this report we propose two al- 
gorithms to solve the distributed slicing problem. The first algorithm improves upon an existing 
algorithm that is based on gossip-based sorting of a set of uniform random numbers. We speed up 
convergence via a heuristic for gossip peer selection. The second algorithm is based on a different 
approach: statistical approximation of the rank of nodes in the ordering. The scalability, efficiency 
and resilience to dynamics of both algorithms relies on their gossip-based models. We present theo- 
retical and experimental results to prove the viability of these algorithms. 

Key-words: Slicing, Gossip, Slice, Churn, Peer-to-Peer, Aggregation, Large Scale, Resource Allo- 
cation. 
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Morcellement distribue dans les systemes dynamiques 



Resume : Un service de morcellement d'un reseau pair-a-pair permet de partitionner les noeuds 
du systeme en plusieurs groupes appeles morceaux. Ce rapport presente deux algorithmes pour 
resoudre le probleme du morcellement reparti. Le premier algorithme ameliore un algorithme exis- 
tant en accerant son temps de convergence. Le second algorithme utilise une approche differente 
d'appoximation statistique. Des resultats theoriques et experimentaux montrent la viabUite de nos 
algorithmes. 

Mots-cles : Morcellement, Bavardage, Morceau, Va-et-vient, Pair-a-pair, Aggregation, Grande 
echelle, Allocation de ressources. 
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1 Introduction 

1.1 Context and Motivations 

The peer to peer (P2P) communication paradigm has now become the prevalent model to build large- 
scale distributed applications, able to cope with both scalability and system dynamics. This is now a 
mature technology: peer to peer systems are slowly moving from application-specific architectures 
to a generic-service oriented design philosophy. More specifically, peer to peer protocols hold the 
promise to integrate into platforms on top of which several applications with various requirements 
may cohabit. This leads to the interesting issue of resource assignment or how to allocate a set of 
nodes for a given application. Examples of targeted platforms for such a service are telecommunica- 
tion platforms, where some set of peers may be automatically assigned to a specific task depending 
on their capabilities, testbed platform such as Planetlab |2|, or desktop-grid-like applications IH. 

In this context, the ordered slicing service has been recently proposed as a building block to 
allocate resources, i.e a set of nodes sharing some characteristics with respect to a given metric 
or attribute, in a large-scale peer to peer system. This service acknowledges the fact that peers 
potentially offer heterogeneous capabilities as revealed by many recent works describing heavy- 
tailed distributions of storage space, bandwidth, and uptime of peers lfT6l [3l [TtI . The slicing service 
ifTOl enables peers in a large-scale unstructured network to self-organize into a partitioning, where 
partitions (slices) are connected overlay networks that represent a given percentage of some resource. 
Such slices can be allocated to specific applications later on. The slicing is ordered in the sense that 
peers get ranked according to their capabilities expressed by an attribute value. 

Large scale dynamic distributed systems consist of many participants that can join and leave at 
will. Identifying peers in such systems that have a similar level of power or capability (for instance, 
in terms of bandwidth, processing power, storage space, or uptime) in a completely decentralized 
manner is a difficult task. It is even harder to maintain this information in the presence of churn. Due 

*Umversidad Rey Juan Carlos, 28933 Mostoles, Spain, anto@gsyc.escet.urjc.es 

tiRISA, INRIA Universite Rennes 1 (ASAP Research Group) 35042 Rennes, France. 
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"'"Uiuversidad Politecnica de Madrid, 28031 Madrid, Spain, emes@eui.upm.es 
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to the intrinsic dynamics of contemporary peer to peer systems it is impossible to obtain accurate 
information about the capabiUties (or even the identity) of the system participants. Consequently, 
no node is able to maintain accurate information about all the nodes. This disqualifies centralized 
approaches. 

Taking this into account, we can summarize the ordered slicing problem we tackle in this report: 
we need to rank nodes depending on their capability, slice the network depending on these capabili- 
ties and, most importantly, readapting the slices continuously to cope with system dynamism. 

Building upon the work on ordered distributed slicing proposed in fTOl, here we focus on the 
issue of accurate slicing. That is, we focus on improving the quality and stability of the slices, both 
aspects being crucial for potential appUcations. 

1.2 Contributions 

The report presents two gossip-based solutions to slice the nodes according to their capability (re- 
flected by an attribute value) in a distributed manner with high probability. The first contribution 
of the report builds upon the distributed slicing algorithm proposed in |10| that we call the JK al- 
gorithm in the sequel of this report. The second algorithm is a different approach based on rank 
approximation through statistical sampling. 

In JK, each node i maintains a random number r^, picked up uniformly at random (between 
and 1), and an attribute value a^, expressing its capability according to a given metric. Each peer 
periodically gossips with another peer j, randomly chosen among the peers it knows about. If the 
order between Vj and is different than the order between aj and a^, random values are swapped 
between nodes. The algorithm ensures that eventually the order on the random values matches the 
order of the attribute ones. The quality of the ranking can then be measured by using a global 
disorder measure expressing the difference between the exact rank and the actual rank of each peer 
along the attribute value. 

The first contribution of this report is to propose a local disorder measure so that a peer chooses 
the neighbor to communicate with in order to maximize the chance of decreasing the global disorder 
measure. The interest of this approach is to speed the convergence up. We provide the analysis and 
experimental results of this improvement. 

Once peers are ordered along the attribute values, the slicing in JK takes place as follows. Ran- 
dom values are used to calculate which slice a node belongs to. For example, a slice containing 
20% of the best nodes according to a given attribute, will be composed of the nodes that end up 
holding random values greater than 0.8. The accuracy of the slicing (independent from the accuracy 
of the ranking) fully depends on the uniformity of the random value spread between and 1 and the 
fact that the proportion of random values between 0.8 and 1 is approximately (but usually not ex- 
actly) 20% of the nodes. Another contribution of this report is to precisely characterize the potential 
inaccuracy resulting from this effect. 

This observation means that the problem of ordering nodes based on uniform random values is 
not fully sufficient for determining slices. This motivates us to find an alternative approach to this 
algorithm and JK in order to determine more precisely the slice each node belongs to. 
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Another motivation for an alternative approach is related to churn and dynamism. It may well 
happen that the churn is actually correlated to the attribute value. For example, if the peers are 
sorted according to their connectivity potential, a portion of the attribute space (and therefore the 
random value space) might be suddenly affected. New nodes will then pick up new random values 
and eventually the distribution of random values will be skewed towards high values. 

The second contribution is an alternative algorithm solving these issues by approximating the 
rank of the nodes in the ordering locally, without the application of random values. The basic idea 
is that each node periodically estimates its rank along the attribute axis depending of the attributes it 
has seen so far. This algorithm is robust and lightweight due to its gossip-based communication pat- 
tern: each node communicates periodically with a restricted dynamic neighborhood that guarantees 
connectivity and provides a continuous stream of new samples. Based on continuously aggregated 
information, the node can determine the slice it belongs to with a decreasing error margin. We show 
that this algorithm provides accurate estimation at the price of a slower convergence. 

1.3 Outline 

The rest of the report is organized as follows: Section |2] surveys some related work. The system 
model is presented in Section |3] The first contribution of an improved ordered slicing algorithm 
based on random values is presented in Section|4]and the second algorithm based on dynamic ranking 
in Section |5] Section |6]concludes the report. 

2 Related Work 

Most of the solutions proposed so far for ordering nodes come from the context of databases, where 
parallelizing query executions is used to improve efficiency. A large majority of the solutions in this 
area rely on centralized gathering or all-to-all exchange, which makes them unsuitable for large-scale 
networks. For instance, the external sorting problem |5 | consists in providing a distributed sorting 
algorithm where the memory space of each processor does not necessarily depend on the input. 
This algorithm must output a sorted sequence of values distributed among processors. The solution 
proposed in [31 needs a global merge of the whole information, and thus it implies a centralization 
of information. Similarly, the percentile finding problem |8|, which aims at dividing a set of values 
into equally sized sets, requires a logarithmic number of all-to-all message exchanges. 

Other related problems are the selection problem and the i/i-quantile search. The selection prob- 
lem |6 , 4J aims at determining the i*'' smallest element with as few comparisons as possible. The 
(j)-quantile search (with e (0, 1]) is the problem to find among n elements the {(pnY'^ element. 
Even though these problems look similar to our problem, they aim at finding a specific node among 
all, while the distributed slicing problem aims at solving a global problem where each node main- 
tains a piece of information. Additionally, solutions to the quantile search problem like the one 
presented in fT3l use an approximation of the system size. The same holds for the algorithm in fT5l, 
which uses similar ideas to determine the distribution of a utility in order to isolate peers with high 
capability — i.e., super-peers. 
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As far as we know, the distributed slicing problem was studied in a P2P system for the first time 
in JTOl. In this report, a node with the fc*'' smallest attribute value, among those in a system of size 
n, tries to estimate its normalized index k/n. The JK algorithm proposed in [10] works as follows. 
Initially, each node draws independently and uniformly a random value in the interval (0, 1] which 
serves as its first estimate of its normalized index. Then, the nodes use a variant of Newscast fVT\ to 
gossip among each other to exchange random values when they find that the relative order of their 
random values and that of their attribute values do not match. This algorithm is robust in face of 
frequent dynamics and guarantees a fast convergence to the same sequence of peers with respect to 
the random and the attribute values. At every point in time the current random value of a node serves 
to estimate the slice to which it belongs (its slice). 

3 Model 

3.1 System model 

We consider a system S containing a set of n uniquely identified nodes. (Value n may vary over 
time, dynamics is explained below). The set of identifiers is denoted by /. Each node can leave 
and new nodes can join the system at any time, thus the number of nodes is a function of time. 
Nodes may also crash. In this report, we do not differentiate between a crash and a voluntary node 
departure. 

Each node i maintains an attribute value a^, reflecting the node capability according to a specific 
metric. These attribute values over the network might have an arbitrary skewed distribution. Initially, 
a node has no global information neither about the structure or size of the system nor about the 
attribute values of the other nodes. 

We can define a total ordering over the nodes based on their attribute value, with the node iden- 
tifier used to break ties. Formally, we let i precede j if and only if ai < aj, or = aj and i < j. 
We refer to this totally ordered sequence as the attribute-based sequence, denoted by A. sequence. 
The attribute-based rank of a node i, denoted hy a.^ e {1, n}, is defined as the index of in 
A. sequence. For instance, let us consider three nodes: 1, 2, and 3, with three different attribute 
values fli = 50, a2 = 120, and as = 25. In this case, the attribute-based rank of node 1 would be 
ai — 2. In the rest of the report, we assume that nodes are sorted according to a single attribute and 
that each node belongs to a unique slice. The sorting along several attributes is out of the scope of 
this report. 

3.2 Distributed Slicing Problem 

Let 5; ti denote the slice containing every node i whose normalized rank, namely — , satisfies I < 
^ < u where I € [0, 1) is the slice lower boundary and u G (0, 1] is the slice upper boundary so that 
all slices represent adjacent intervals {li, ui], {I2, U2]... Let us assume that we partition the interval 
(0, 1] using a set of slices, and this partitioning is known by aU nodes. The distributed slicing 
problem requires each node to determine the slice it currently belongs to. Note that the problem 
stated this way is similar to the ordering problem, where each node has to determine its own index 
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Figure 1: Slicing of a population based on a height attribute. 



in A. sequence. However, the reference to slices introduces special requirements related to stability 
and fault tolerance, besides, it allows for future generalizations when one considers different types 
of categorizations. 

Figure [T] illustrates an example of a population of 10 persons, to be sorted against their height. 
A partition of this population could be defined by two slices of the same size: the group of short 
persons, and the group of tall persons. This is clearly an example where the distribution of attribute 
values is skewed towards 2 meters. The rank of each person in the population and the two slices are 
represented on the bottom axis. Each person is represented as a small cross on these axesQ Each 
slice is represented as an oval. The slice 5*1 = S^ i contains the five shortest persons and the sUce 
S2 = Si I contains the five tallest persons. 

Observe that another way of partitioning the population could be to define the group of short 
persons as that containing all the persons shorter than a predefined measure (e.g., 1.65m) and the 
group of tall persons as that containing the persons taller than this measure. However, this way of 
partitioning would most certainly lead to an unbalanced distribution of persons, in which, for in- 
stance a group might be empty (while a slice is almost surely non-empty). Since the distribution of 
attribute values is unknown and hard to predict, defining relevant groups is a difficult task. For exam- 
ple, if the distribution of the human heights were unknown, then the persons taller than Im could be 
considered as tall and the persons shorter than Im could be considered as short. Conversely, slices 
partition the population into subsets representing a predefined portion of this population. Therefore, 
in the rest of the report, we consider slices as defined as a proportion of the network. 

3.3 Facing Churn 

Node churn, that is, the continuous arrival and departure of nodes is an intrinsic characteristic of peer 
to peer systems and may significantly impact the outcome, and more specifically the accuracy of the 
slicing algorithm. The easier case is when the distribution of the attribute values of the departing 
and arriving nodes are identical. In this case, in principle, the arriving nodes must find their slices, 
but the nodes that stay in the system are mostly able to keep their slice assignment. Even in this 
case however, nodes that are close to the border of a slice may expect frequent changes in their slice 
due to the variance of the attribute values, which is non-zero for any non-constant distribution. If 
the arriving and departing nodes have different attribute distributions, so that the distribution in the 
actual network of live nodes keeps changing, then this effect is amplified. However, we believe that 

'Note that the shortest (resp. largest) rank is represented by a cross at the extreme left (resp. right) of the bottom axis. 
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this is a realistic assumption to consider that the churn may be correlated to some specific values (for 
example if the considered attribute is uptime or connectivity). 

4 Dynamic Ordering by Exchange of Random Values 

This section proposes an algorithm for the distributed slicing problem improving upon the original 
JK algorithm [ lOJ, by considering a local measure of the global disorder function. In this section we 
present the algorithm along with the corresponding analysis and simulation results. 

4.1 On Using Random Numbers to Sort Nodes 

This Section presents the algorithm built upon JK. We refer to this algorithm as mod-JK (standing 
for modified JK). In JK, each node i generates a number ri E (0, 1] independently and uniformly 
at random. The key idea is to sort these random numbers with respect to the attribute values by 
swapping these random numbers between nodes, so that if < aj then < Vj . Eventually, the 
attribute values (that are fixed) and the random values (that are exchanged) should be sorted in the 
same order. That is, each node would like to obtain the a;*'* largest random number if it owns the x*^ 
largest attribute value. Let R. sequence denote the random sequence obtained by ordering all nodes 
according to their random number. Let pi{t) denote the index of node i in R. sequence at time t. 
When not required, the time parameter is omitted. 

To illustrate the above ideas, consider that nodes 1, 2, and 3 from the previous example have 
three distinct random values: ri = 0.85, r2 — 0.1, and — 0.35. In this case, the index pi of node 
1 would be 3. Since the attribute values are — 50, 02 — 120, and 03 = 25, the algorithm must 
achieve the following final assignment of random numbers: ri — 0.35, r2 = 0.85, and = 0.1. 

Once sorted, the random values are used to determine the portion of the network a peer belongs 

to. 

4.2 Definitions 

View. Every node i keeps track of some neighbors and their age. The age of neighbor j is a 
timestamp, tj, set to when j becomes a neighbor of i. Thus, node i maintains an array containing 
the id, the age, the attribute value, and the random value of its neighbors. This array, denoted Mi, is 
called the view of node i. The views of all nodes have the same size, denoted by c. 

Misplacement. A node participates in the algorithm by exchanging its rank with a misplaced 
neighbor in its view. Neighbor j is misplaced if and only if 

• ai > aj and ri < rj , or 

• Qi < aj and 7\ > Vj . 

We can characterize these two cases by the predicate {aj — a;) (r^ — r^) < 0. 
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Global Disorder Measure. In ifTOl . a measure of the relative disorder of sequence R. sequence 
with respect to sequence A. sequence was introduced, called the global disorder measure (GDM) 
and defined, for any time t, as 

GDM{t) = iV(«^-p(t).)'. 
n ^ — ' 

i 

The minimal value of GDM is 0, which is obtained when p{t)i = OLi for all nodes i. In this case 
the attribute-based index of a node is equal to its random value index, indicating that random values 
are ordered. 



4.3 Improved Ordering Algorithm 

In this algorithm, each node i searches its own view Mi for misplaced neighbors. Then, one of them 
is chosen to swap random value with. This process is repeated until there is no global disorder In 
this version of the algorithm, we provide each node with the capability of measuring locally the 
disorder. This leads to a new heuristic for each node to determine the neighbor to exchange with 
which decreases most the disorder 

The proposed technique attempts to decrease the global disorder in each exchange as much as 
possible via selecting the neighbor from the view that minimizes the local disorder (or, equivalently, 
maximizes the order gain) as defined below. 

For a node i to evaluate the gain of exchanging with a node j of its current view Mi, we define 
its local disorder measure (abbreviated LDMi). Let LA.sequence^ and LR.sequence^ be the local 
attribute sequence and the local random sequence of node i, respectively. These sequences are 
computed locally by i using the information A/i U {i}. Similarly to A. sequence and R. sequence, 
these are the sequences of neighbors where each node is ordered according to its attribute value and 
random number, respectively. Let, for any j E MiU {i}, £pj{t) and £aj{t) be the indices of rj and 
Qj in sequences LR.sequence^^ and LA. sequence. i, respectively, at time {t). At any time t, the local 
disorder measure of node i is defined as: 

LDM,{t) = -1^ (^"^■(*) - ^Z'^- 

je7V.(t)u{i} 

We denote by Gij{t + 1) the reduction on this measure that i obtains after exchanging its random 
value with node j between time t and t + 1. We define it as: 

G,j(t +1) = LDMdt) ~ LDM,{t + 1), 

c+1 

The heuristic used chooses for node i the misplaced neighbor j that maximizes Gi j + 1). 
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Variable 


Description 


3 


the identifier of tlie neigiibor 


ij 


tlie age of the neighbor 


aj 


the attribute value of the neighbor 




the random value of the neighbor 



Table 1: The array corresponding to the view entry of the neighbor (j £ J\fi). 



4.3.1 Sampling Uniformly at Random 

The algorithm relies on the fact that potential misplaced nodes are found so that they can swap 
their random numbers thereby increasing order. If the global disorder is high, it is very likely that 
any given node has misplaced neighbors in its view to exchange with. Nevertheless, as the system 
gets ordered, it becomes more unlikely for a node i to have misplaced neighbors. In this stage the 
way the view is composed plays a crucial role: if fresh samples from the network are not available, 
convergence can be slower than optimal. 

Several protocols may be used to provide a random and dynamic sampling in a peer to peer 
system such as Newscast |12|, Cyclon ||18 I or Lpbcast [9J. They differ mainly by their closeness to 
the uniform random sampling of the neighbors and the way they handle churn. In this report, we 
chose to use a variant of the Cyclon protocol to construct and update the views Q, as it is reportedly 
the best approach to achieve a uniform random neighbor set for all nodes. 

4.3.2 Description of the Algorithm 

The algorithm is presented in Figure |2] The active thread at node i runs the membership (gossip- 
ing) procedure (recompute-view()i) and the exchange of random values periodically. As motivated 
above, the membership procedure, specified in Figure [3] is similar to the Cyclon algorithm: each 
node i maintains a view A/i containing one entry per neighbor. The entry of a neighbor j corre- 
sponds to a tuple presented in Table [1] Node i copies its view, selects the oldest neighbor j of its 
view, removes the entry Cj of j from the copy of its view, and finally sends the resulting copy to 
j. When j receives the view, j sends its own view back to i discarding possible pointers to i, and i 
and j update their view with the one they receive. This variant of Cyclon, as opposed to the original 
version, exchanges all entries of the view at each step. 

The algorithm for exchanging random values from node i starts by measuring the ordering that 
can be gained by swapping with each neighbor (Lines |4]-[8]l. Then, i chooses the neighbor j S A/i 
that maximizes gain ^ for any of its neighbor k. Formally, i finds j S A/i such that for any 
k £ JVi,we have 

+ > G,,fe(t + 1). (2) 
Using the definition of dj in Equation ([l]). Equation (|2|i is equivalent to 

£a^{t)£pj{t) + eaj{t)ip,{t) ~ £aj{t)epj{t) > iar{t)£pk{t) + £ak{t)£p,{t) ~ £ak{t)£pk{t). 



INRIA 



Distributed Slicing in Dynamic Systems 



11 



Initial state of node i 

(1) period ^, initially set to a constant; 

Ti, arandom value chosen in (0, 1]; a^, the attribute value; 

slicei <— J_, the slice i belongs to; A/j, the view; 

gaiUji, a real value indicating the gain achieved by exchanging with j'; 

gain-max = 0, a real. 

Active thread at node i 

(2) wa\t{periodi) 

(3) recompute-view()j 

(4) forj'eM 

(5) if gaiUji > gain-max then 

(6) gain-max ^ gainy 

(7) j «- f 

(8) end for 

(9) send(REQ,rj,a,) to j 

(10) recv(ACK,rp from j 

(11) rj^rr 

(12) if {aj - ai)(rj - r^) < then 

(13) n ^ rj 

(14) slicei such that I < ri < u 

Passive thread at node i activated upon reception 

(15) recv(REQ, rj,aj) from j 

(16) send(ACK,ri)toj 

(17) if (aj — ai){rj — r,) < then 

(18) n Tj 

(19) slicei <— <5j,u such that I < ri <u 



Figure 2: Dynamic ordering by exchange of random values. 
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Active thread at node i 

(1) for j' e Ni do tj, ^ tj, + 1 end for 

(2) 3 ^ j" ■ tj" = m3Xj'g^.(tj/) 

(3) send(REQ',A/'i \ {ej} U {(i, 0, a,, r,)}) to j 

(4) recv{ACK',A^j ) from j 

(5) duplicated- entries = {e : e.id G J\fj nA/";} 

(6) Afi <— Afi U (A/j \ duplicated- entries \ {ci}) 

Passive tliread at node i activated upon reception 

(7) recv {REQ',Afj) from j 

(8) send(ACK',M)toj 

(9) duplicated- entries = {e £ Afj : e.id £ Afj nAfi} 

(10) Afi <— Afi U (A/j \ duplicated- entries) 



Figure 3: recompute-view(): procedure used to update the view based on a simple variant of the 
Cyclon algorithm. 

In Figure|2]of node i, we refer to gaiUj as the value of £ai{t)£pj{t) + £aj{t)£pi{t) — £aj{t)£pj{t). 

From this point on, i exchanges its random value with the random value rj of node j (LinefTTTi. 
The passive threads are executed upon reception of a message. In Figure |2] when j receives the 
random value of node i, it sends back its own random value Vj for the exchange to occur (LinesfTSl- 
[T6b . Observe that the attribute value of i is also sent to j, so that j can check if it is correct to 
exchange before updating its own random number (Lines [TTl - frsl l. Node i does not need to receive 
attribute value aj of j, since i already has this information in its view and the attribute value of a 
node never changes over time. 

4.4 Analysis of Slice Misplacement 

In mod-JK, as in JK, the current random number of a node i determines the slice Si of the node. 
The objective of both algorithms is to reduce the global disorder as quickly as possible. Algorithm 
mod-JK consists in choosing one neighbor among the possible neighbors that would have been 
chosen in JK, plus the GDM of JK has been shown to fit an exponential decrease ifTOl . Consequently 
mod-JK experiences also an exponential decrease of the global disorder Eventually, JK and mod-JK 
ensure that the disorder has fully disappeared. However, the accuracy of the slices heavily depends 
on the uniformity of the random value spread between and 1 . It may happen, that the distribution 
of the random values is such that some peers decide upon a wrong slice. Even more problematic is 
the fact that this situation is unrecoverable unless a new random value is drawn for all nodes. This 
may be considered as an inherent limitation of the approach. For example, consider a system of size 
2, where nodes 1 and 2 have the random values ri = 0.1, r2 = 0.4. If we are interested in creating 
two slices of equal size, the first slice will be of size 2 and the second of size zero, even after perfect 
ordering of the random values. 

Therefore, an important step is to characterize the inaccuracy of the uniform distribution to access 
the potential impact on the slice assignment resulting from the fact that uniformly generated random 
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numbers are not distributed perfectly evenly throughout the domain. First of all, consider a slice Sp 
of length p. In a network of n nodes, the number of nodes that will fall into this slice is a random 
variable X with a binomial distribution with parameters n and p. The standard deviation of X is 
therefore A/rip(T^^p). This means that the relative proportional expected difference from the mean 
(i.e., np) can be approximated as ^(1 — p)/{np), which is very large if p is small, in fact, goes to 
infinity as p tends to zero, although a very large n compensates for this effect. For a "normal" value 
of p, and a reasonably large network, the variance is very low however. 

To stay with this random variable, the following result bounds, with high probability, its deviation 
from its mean. 

Lemma 4.1. For any [3 € (0, 1], a slice Sp of length p £ (0, 1] has a number of peers X G 
[(1 — (3)np, (1 + I3)np] with probability at least 1 — e as long as p > ln(2/e). 

Proof. The way nodes choose their random number is like drawing n times, with replacement and 
independently uniformly at random, a value in the interval (0, 1]. Let Xi, ...,X„ be the n corre- 
sponding independent identically distributed random variables such that: 

Xi = 1 if the value drawn by node i belongs to Sp and 
Xi ~ otherwise. 

We denote X — J^^^i number of elements of interval Sp drawn among the n drawings. 

The expectation of X is np. From now on we compute the probability that a bounded portion of the 
expected elements are misplaced. Two Chernoff bounds [14] give: 



Pi[\X - np\ > I3np] < 2e~ 



Pr[X > {1 + I3)np] <e-^ \ 
Pt[X < (1 - /3)?ip] < J 
with < /3 < 1. That is, the probability that more than (/3 time the number expected) elements are 



misplaced regarding to interval Sp is bounded by 2e 3 . We want this to be at most e. This yields 
the result. □ 

To measure the effect discussed above during the simulation experiments, we introduce the slice 
disorder measure (SDM) as the sum over all nodes i of the distance between the slice i actually 
belongs to and the slice i believes it belongs to. For example (in the case where all slices have the 
same size), if node i belongs to the slice (according to its attribute value) while it thinks it belongs 
to the 3'^'' slice (according to its rank estimate) then the distance for node iis |1 — 3| — 2. Formally, 
for any node i, let S^.u be the actual correct slice of node i and let S^, f, (t) be the slice i estimates 
as its slice at time t. The slice disorder measure is defined as: 



SDMit)^J2 — 

^ ^ ?/,- — 



Ui li 

i 

SDM{t) is minimal (equals 0) if for all nodes i, we have S^ f (t) = S^^ii. 

In fact, it is simple to show that, in general, the probability of dividing n peers into two slices of 
the same size is less than y^2/mT. This value is very small even for moderate values of n. Hence, it 
is highly possible that the random number distribution does not lead to a perfect division into slices. 
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4.5 Simulation Results 

We present simulation results using PeerSim ifTTl . using a simplified cycle-based simulation model, 
where all messages exchanges are atomic, so messages never overlap. First, we compare the perfor- 
mance of the two algorithms: JK and mod-JK. Second, we study the impact of concurrency that is 
ignored by the cycle-based simulations. 

4.5.1 Performance Comparison 

We compare the time taken by these algorithms to sort the random values according to the attribute 
values (i.e., the node with the largest attribute value of the system value obtains the random 
value). In order to evaluate the convergence speed of each algorithm, we use the slice disorder 
measure as defined in Section l4~4l 

We simulated 10^ participants in 100 equally sized slices (when unspecified), each with a view 
size c — 20. Figure [4(a)] illustrates the difference between the global disorder measure and the slice 
disorder measure while Figure [4(b)| presents the evolution of the slice disorder measure over time 
for JK, and mod-JK. 

Figure [4(a)] shows the different speed at which the global disorder measure and the slice disorder 
measure converge. When values are sufficiently large, the GDM and SDM seem tightly related: if 
GDM increases then SDM increases too. Conversely, there is a significant difference between the 
GDM and SDM when the values are relatively low: the GDM reaches while the SDM is lower 
bounded by a positive value. This is because the algorithm does lead to a totally ordered set of 
nodes, while it still does not associate each node with its correct slice. Consequently the GDM is 
not sufficient to rightly estimate the performance of our algorithms. 

Figure [4(b)| shows the slice disorder measure to compare the convergence speed of our algorithm 
to that of JK with 10 equally sized slices. Our algorithm converges significantly faster than JK. Note 
that none of the algorithm reaches zero SDM, since they are both based on the same idea of sorting 
randomly generated values. Besides, since they both used an identical set of randomly generated 
values, both converge to the same SDM. 

4.5.2 Concurrency 

The simulations are cycle-based and at each cycle an algorithm step is done atomically so that 
no other execution is concurrent. More precisely, the algorithms are simulated such that in each 
cycle, each node updates its view before sending its random value or its attribute value. Given this 
implementation, the cycle-based simulator does not allow us to realistically simulate concurrency, 
and a drawback is that view is up-to-date when a message is sent. In the following we artificially 
introduce concurrency (so that view might be out-of-date) into the simulator and show that it has 
only a slight impact on the convergence speed. 

Introducing concurrency might result in some problems because of the potential staleness of 
views: unsuccessful swaps due to useless messages. Technically, the view of node i might indicate 
that i has a random value r while this value is no longer up-to-date. This happens if i has lastly 
updated its view before j swapped its random value with another f . Moreover, due to asynchrony. 
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Figure 4: (a) Evolution of the global disorder measure over time, (b) Slice disorder measure over 
time, (c) Percentage of unsuccessful swaps in the ordering algorithms, (d) Convergence speed under 
high concurrency. 
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it could happen that by the time a message is received this message has become useless. Assume 
that node i sends its random value Vi to j in order to obtain Vj at time t and j receives it by time 
t + 5. With no loss of generality assume > rj. Then if j swaps its random value with / such that 
r-j > r, between time t and t + 5, then the message of i becomes useless and the expected swap does 
not occur (we call this an unsuccessful swap). 

Figure 4(d)| indicates the impact of concurrent message exchange on the convergence speed while 
Figure [4(c) shows the amount of useless messages that are sent. Now, we explain how the concur- 
rency is simulated. Let the overlapping messages be a set of messages that mutually overlap: it 
exists, for any couple of overlapping messages, at least one instant at which they are both in-transit. 
For each algorithm we simulated ( i) full concurrency: in a given cycle, all messages are overlapping 
messages; and (ii) half concurrency: in a given cycle, each message is an overlapping message with 
probability i. Generally, we see that increasing the concurrency increases the number of useless 
messages. Moreover, in the modified version of JK, more messages are ignored than in the original 
JK algorithm. This is due to the fact that some nodes (the most misplaced ones) are more likely tar- 
geted which increases the number of concurrent messages arriving at the same nodes. Since a node i 
ignored more likely a message when it receives more messages during the same cycle, it comes out 
that concentrating message sending at some targets increases the number of useless messages. 

Figure |4(d)| compares the convergence speed under full concurrency and no concurrency. We 
omit the curve of half-concurrency since it would have been similar to the two other curves. Full- 
concurrency impacts on the convergence speed very slightly. 



5 Dynamic Ranking by Sampling of Attribute Values 

In this section we propose an alternative algorithm for the distributed slicing problem. This algorithm 
circumvents some of the problems identified in the previous approach by continuously ranking nodes 
based on observing attribute value information. Random values no longer play a role, so non-perfect 
uniformity in the random value distribution is no longer a problem. Besides, this algorithm is not 
sensitive to churn even if it is correlated with attribute values. 

In the remaining part of the report we refer to this new algorithm as the ranking algorithm while 
referring to JK and mod-JK as the ordering algorithms. Here, we elaborate on the drawbacks arising 
from the ordering algorithms relying on the use of random values that are solved by the ranking 
approach. 

Impact of attribute correlated with dynamics. As already mentioned, the ordering algorithms 
rely on the fact that random values are uniformly distributed. However, if the attribute values are 
not constant but correlated with the dynamic behavior of the system, the distribution of random 
values may change from uniform to skewed quickly. For instance, assume that each node maintains 
an attribute value that represents its own lifetime. Although the algorithm is able to quickly sort 
random values, so nodes with small lifetime will obtain the small random values, it is more likely 
that these nodes leave the system sooner than other nodes. This results in a higher concentration of 
high random values and a large population of the nodes wrongly estimate themselves as being part 
of the higher slices. 
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Inaccurate slice assignments. As discussed in previous sections in detail, slice assignments will 
typically be imperfect even when the random values are perfectly ordered. Since the ranking ap- 
proach does not rely on ordering random nodes, this problem is not raised: the algorithm guarantees 
eventually perfect assignment in a static environment. 

Concurrency side-effect. In the previous ordering algorithms, a non negligible amount of mes- 
sages are sent unnecessarily. The concurrency of messages has a drastic effect on the number of 
useless messages as shown previously, slowing down convergence. In the ranking algorithm concur- 
rency has no impact on convergence speed because all received messages are taken in account. This 
is because the information encapsulated in a message (the attribute value of a node) is guaranteed to 
be up to date, as long as the attribute values are constant, or at least change slowly. 

5.1 Ranking Algorithm Specification 

The pseudocode of the ranking algorithm is presented in Figure |5] As opposed to the ordering 
algorithm of the previous section, the ranking algorithm does not assign random initial unalterable 
values as candidate ranks. Instead, the ranking algorithm improves its rank estimate each time a new 
message is received. 

The ranking algorithm works as follows. Periodically each node i updates its view Mi following 
an underlying protocol that provides a uniform random sample (Line [3]i; later, we simulate the 
algorithm using the variant of Cyclon protocol presented in Section l4.3.2l Node i computes its rank 
estimate (and hence its slice) by comparing the attribute value of its neighbors to its own attribute 
value. This estimate is set to the ratio of the number of nodes with a lower attribute value that i 
has seen over the total number of nodes i has seen (Line[T5]l. Node i looks at the normalized rank 
estimate of all its neighbors. Then, i selects the node ji closest to a slice boundary (according to the 
rank estimates of its neighbors). Node i selects also a random neighbor j2 among its view (Line[T2]l. 
When those two nodes are selected, i sends an update message, denoted by a flag UPD, to ji and j2 
containing its attribute value (Line[T3}{T4|l. 

The reason why a node close to the slice boundary is selected as one of the contacts is that such 
nodes need more samples to accurately determine which slice they belong to (subsection l5.2l shows 
this point). This technique introduces a bias towards them, so they receive more messages. 

Upon reception of a message from node i, the passive threads of ji and j2 are activated so that 
ji and i2 compute their new rank estimate ^ and rj^ . The estimate of the slice a node belongs to, 
foUows the computation of the rank estimate. Messages are not replied, communication is one-way, 
resulting in identical message complexity to JK and mod-JK. 

5.2 Theoretical Analysis 

The following Theorem shows a lower bound on the probability for a node i to accurately estimate 
the slice it belongs to. This probabiUty depends not only on the number of attribute exchanges but 
also on the rank estimate of i. 
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Initial state of node i 

(1) period^, initially set to a constant; rj, a value in (0, 1] ; 
ai, the attribute value; 6, the closest slice boundary to node i; 
Qi, the counter of encountered attribute values; li, the counter 
of lower attribute values; slicei *— ±; A/i, the view. 

Active thread at node i 

(2) wait(perio(ii ) 

(3) recompute-view()i 

(4) dist-min ^ oo 

(5) forj'eM 

(6) Si ^ Si + 1 

(7) if aj/ < ai then £i ^ £{ -\- 1 

(8) if dist(aj/ , 6) < dist-min then 

(9) dist-min <— dist(aj/ , &) 

(10) ji ^ j' 

(11) end for 

(12) Let j2 be a random node of A/i 

(13) send(UPD, ai) to ji 

(14) send{UPD,ai) toj2 

(15) n ^ ii/gi 

(16) slice <— Si^y, such that I < ri < u 

Passive thread at node i activated upon reception 

(17) recv(UPD,aj) from j 

(18) if aj < ai then^i <- 4 H- 1 

(19) gi^9^ + l 

(20) Vi ^ e,/g, 

(21) slice ^ such that I < < u 



Figure 5: Dynamic ranking by exchange of attribute values. 
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Theorem 5.1. Let p be the normalized rank of i and let p be its estimate. For node i to exactly 
estimate its slice with confidence coefficient of 100{l — 0;)%, the number of messages i must receive 
is: 



where d is the distance between the rank estimate of i and the closest slice boundary, and 
represents the endpoints of the confidence interval. 

Proof. Each time a node receives a message, it checks whether or not the attribute value is larger 
or lower than its own. Let Xi, ...,Xk be fc (fc > 0) independent identically distributed random 
variables described as follows. Xj = 1 with probability ^ = P (indicating that the attribute value is 
lower) and j e {1, fc}, otherwise Xj = (indicating the attribute value is larger). By the central 
Umit theorem, we assume fc > 30 and we approximate the distribution of X ^ Sj=i -^j '^he 

normal distribution. We estimate X hy X — X]j=i -^j ^'^'^ phy p = ^. 

We want a confidence coefficient with value 1 — a. Let $ be the standard normal distribution 
function, and let be <i>^^(l — ^). Now, by the Wald large-sample normal test in the binomial 



case, where the standard deviation of p is a^p) — ^^'^ , we have: 



p-p 



< Z2 



a{p) 

p~Z^a{p) <p <p + Z!^a{p). 
Next, assume that p falls into the slice S'; „, with / and u its lower and upper boundaries, respec- 



tively. Then, as long as p — Z^ \J ^^^j^ > I and p + Z^ y ^'^^^ < u, the slice estimate is exact 
with a confidence coefficient of 100(1 — a)%. Let d ~ min(p — l,u — p), then we need 



> 



□ 

To conclude, under reasonable assumptions all node estimate its slice with confidence coefficient 
100(1 — a)%, after a finite number of message receipts. Moreover a node closer to the slice boundary 
needs more messages than a node far from the boundary. 



5.3 Simulation Results 

This section evaluates the ranking algorithm by focusing on three different aspects. First, the per- 
formance of the ranking algorithm is compared to the performance of the ordering algorithrr0 in a 

^We omit comparison witli JK since tlie performance obtained with mod-JK are eitlier similai' or better 
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Figure 6: (a) Comparing performance of the ordering algorithm and the ranking algorithm (static 
case), (b) Comparing the ranking algorithm on top of a uniform drawing or a Cyclon-like protocol. 

(c) Effect of dynamics burst on the convergence of the ordering algorithm and the ranking algorithm. 

(d) Effect of a low and regular churn on the convergence of the ordering algorithm and the ranking 
algorithm. 



large-scale system where the distribution of attribute values does not vary over time. Second, we 
investigate if sufficient uniformity is achievable in reality using a dedicated protocol. Third, the 
ranking algorithm and ordering algorithm are compared in a dynamic system where the distribution 
of attribute values may change. 

For this purpose, we ran two simulations, one for each algorithms. The system contains (initially) 
10^ nodes and each view contains 10 uniformly drawn random nodes and is updated in each cycle. 
The number of slices is 100, and we present the evolution of the slice disorder measure over time. 
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5.3.1 Performance Comparison in the Static Case 

Figure |6(a)| compares the ranking algorithm to the ordering algorithm while the distribution of at- 
tribute values do not change over time (varying distribution is simulated below). 

The difference between the ordering algorithm and the ranking algorithm indicates that the rank- 
ing algorithm gives a more precise result (in terms of node to slice assignments) than the ordering 
algorithm. More importantly, the slice disorder measure obtained by the ordering algorithm is lower 
bounded while the one of the ranking algorithm is not. Consequently, this simulation shows that the 
ordering algorithm might fail in sUcing the system while the ranking algorithm keeps improving its 
accuracy over time. 

5.3.2 Feasibility of the Ranking Algorithm 

Figure [6(b)| shows that the ranking algorithm does not need artificial uniform drawing of neighbors. 
Indeed, an underlying view management protocol might lead to similar performance results. In the 
presented simulation we used an artificial protocol, drawing neighbors randomly at uniform in each 
cycle of the algorithm execution, and the variant of the Cyclon [ 18] view management protocol pre- 
sented above. Those underlying protocols are distinguished on the figure using terms "uniform" (for 
the former one) and "views" (for the later one). As said previously, the Cyclon protocol consists of 
exchanging views between neighbors such that the communication graph produced shares similar- 
ities with a random graph. This figure shows that both cases give very similar results. The SDM 
legend is on the right-handed vertical axis while the left-handed vertical axis indicates what percent- 
age the SDM difference represents over the total SDM value. At any time during the simulation (and 
for both type of algorithms) its value remains within plus or minus 7%. The two SDM curves of the 
ranking algorithm almost overlap. Consequently, the ranking algorithm and the variant of Cyclon 
presented in subection l4.3.2l achieve very similar result. 

To conclude, the variant of Cyclon algorithm presented in the previous section can be used with 
the ranking algorithm to provide the shuffling of views. 

5.3.3 Performance Comparison in the Dynamic Case 

In Figure [6(c)| each of the two curves represents the slice disorder measure obtained over time using 
the ordering algorithm and the ranking algorithm respectively. We simulate the churn such that 
0. 1 % of nodes leave and 0. 1 % of the nodes join in each cycle during the 200 first cycles. We observe 
how the SDM converges. The churn is reasonably and pessimistically tuned compared to recent 
experimental evaluations 1 17| of the session duration in three well-known P2P systems|l 

The distribution of the churn is correlated to the attribute value of the nodes. The leaving nodes 
are the nodes with the lowest attribute values while the entering nodes have higher attribute values 
than all nodes already in the system. The parameter choices are motivated by the need of simulating 
a system in which the attribute value corresponds to the session duration of nodes, for example. 

'in 1171 . roughly all nodes have left the system after 1 day while there are still 50% of nodes after 25 minutes. In our case, 
assuming that in average a cycle lasts one second would lead to more than 54% of leave in 9 minutes. 
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The churn introduces a significant disorder in the system which counters the fast decrease. When, 
the churn stops, the ranking algorithm readapts well the slice assignments: the SDM starts decreasing 
again. However, in the ordering algorithm, the convergence of SDM gets stuck. This leads to a poor 
slice assignment accuracy. 

In Figure |6(d)[ each of the two curves represent the slice disorder measure obtained over time 
using the ordering algorithm, the ranking algorithm, and a modified version of the ranking algorithm 
using attribute values recorded in a sliding-window, respectively. (The simulation obtained using 
sliding windows is described in the next subsection.) The churn is diminished and made more 
regular than in the previous simulation such that 0.1% of nodes leave and 0.1% of nodes join every 
10 cycles. 

The curves fits a fast decrease (superlinear in the number of cycles) at the beginning of the 
simulation. At first cycles, the ordering gain is significant making the impact of churn negligible. 
This phenomenon is due to the fact that SDM decreases rapidly when the system is fully disordered. 
Later on, however, the decrease slope diminishes and the churn effect reduces the amount of nodes 
with a low attribute value while increasing the amount of nodes with a large attribute value. This 
unbalance leads to a messy slice assignment, that is, each node must quickly find its new slice to 
prevent the SDM from increasing. In the ordering algorithm the SDM starts increasing from cycle 
120. Conversely, with the ranking algorithm the SDM starts increasing not earlier than at cycle 730. 
Moreover the increase slope is much larger in the former algorithm than in the latter one. 

Even though the performance of the ranking algorithm are really significant, its adaptiveness to 
churn is not surprising. Unlike the ordering algorithm, the ranking one keeps re-estimating the rank 
of each node depending on the attribute values present in the system. Since the churn increases the 
attribute values present in the system, nodes tend to receive more messages with higher attribute 
values and less messages with lower attribute values, which turns out to keep the SDM low, despite 
churn. Further on, we propose a solution based on sliding-window technique to limit the increase of 
the SDM in the ranking algorithm. 

To conclude, the results show that when the churn is related to the attribute (e.g., attribute rep- 
resents the session duration, uptime of a node), then the ranking algorithm is better suited than the 
ordering algorithm. 

5.3.4 Sliding-window for Limiting the SDM Increase 

In Figure [6(d)] the "sliding-window" curve presents a slightly modified version of the ranking algo- 
rithm that encompasses SDM increase due to churn correlated to attribute values. Here, we present 
this enrichment. 

In Section|5] the ranking algorithm specifies that each node takes into account all received mes- 
sages. More precisely, upon reception of a new message each node i re-computes immediately its 
rank estimate and the slice it thinks it belongs to without remembering the attribute values it has 
seen. Consequently the messages received long-time ago have as much importance as the fresh mes- 
sages in the estimate of i. The drawback, as it appeared in Figure [6(d)| of Section |431 is that if the 
attribute values are correlated to churn, then the precision of the algorithm might diminish. 

To cope with this issue, the previous algorithm can be easily enriched in the following way. 
Upon reception of a message, each node i records an information about the attribute value received 
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in a fixed-size ordered set of values. Say this set is a first-in first-out buffer such that only the most 
recent values remain. Right after having recorded this information, node i can re-compute its rank 
estimate and its slice estimate based on the most relevant piece of information (having discarded the 
irrelevant piece). Consequently, the estimate would rely only on fresh attribute values encountered 
so that the algorithm would be more tolerant to changes (e.g., dynamics or non-uniform evolution of 
attribute values). Of course, since the analysis (cf. Section |S!2] i shows that nodes close to the slice 
boundary require a large number of attribute values for estimating precisely their estimates, it would 
be unaffordable to record all these last attribute values encountered due to space limitation. 

Actually, the only necessary relevant information of a message is simply whether it contains a 
lower attribute value than the attribute value of i, or not. Consequently, a single bit per message 
would be sufficient to record the necessary information (e.g., adding a 1 meaning that the attribute 
value is lower, and otherwise). Thus, even though a node i would require 10'' messages to rightly 
estimate its slice (with high probability), node i simply needs to allocate an array of size 10"*/ (8 * 
1000) = l,25kB. 

As expected. Figure |6(d)| shows that the sliding-window method applied to the ranking algo- 
rithm prevents its SDM from increasing. Consequently, at some point in time, the resulting slice 
assignment may become even more accurate. 

6 Conclusion 
6.1 Summary 

Peer to peer systems may now be turned into general frameworks on top of which several applications 
might cohabit. To this end, allocating resources to applications, according to their needs require 
specific algorithms to partition the network in a relevant way. The ordered slicing algorithm proposed 
in ifTOl provided a first attempt to "slice" the network, taking into account the potential heterogeneity 
of nodes. This algorithm relies on each node drawing a random value uniformly and swapping 
continuously those random values, with candidate nodes, so that the order between attributes values 
(reflecting the capabilities of nodes) and random ones match. Results from flOl have shown that 
slices can be maintained efficiently and in large-scale systems even in the presence of churn. 

In this report, we first proposed an improvement over the initial sorting algorithm based on 
a judicious choice of candidate nodes to swap values. This is based on each node being able to 
estimate locally the potential decrease of the global disorder measure. We provided an analysis 
along with some simulation results showing that the convergence speed is significantly improved. 
We then identified two issues related to the use of static random values. The first one refers to the 
fact that slice assignment heavily depends on the degree of uniformity of the initial random value. 

The second is related to the fact that once sorted along one attribute axis, the churn (or failures) 
might be correlated to the attribute, therefore leading to a unrecoverable skewed distribution of the 
random values resulting in a wrong slice assignment. Our second contribution is an algorithm en- 
abling nodes to continuously re-estimate their rank relatively to other nodes based on their sampling 
of the network. 
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6.2 Perspective 

This report used a variant of the Cyclon protocol to obtain quasi-uniform distribution of neighbors. 
There are various protocols that might be used for different purpose. For instance, Newscast can be 
used for its resilience to very high dynamics as in |10|. Some other protocols exist in the literature. 
Deciding exactly how to parameterize the underlying peer sampling service might be an interesting 
future direction. 
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